Atomic wave packet basis for quantum information 



Ashok Muthukrishnanu and C. R. Stroud, Jr. 
The Institute of Optics, University of Rochester, Rochester, New York 14-627 

(March 15, 2001) 



O 
O 

(N 



> 

in 
O 

o 

22 

9 L,: 
-i— > ■ 

G ■ 
cd 

3 : 
cr 



13 



We propose a wave packet basis for storing and processing several qubits of quantum information in 
a single multilevel atom. Using radially localized wave packet states in the Rydberg atom, we construct 
an orthogonal basis that is related to the usual energy level basis by a quantum Fourier transform. 
A transform-limited laser pulse that is short compared with the classical Kepler period of the system 
interacts mainly with the wave packet state localized near the atomic core, allowing selective control in 
this basis. We argue that wave packet control in this regime is useful for multilevel quantum information 
processing. 

PACS numbers: 03.67.-a, 32.80.Qk 



Quantum information processing is a growing interdis- 
ciplinary field that has seen rapid progress in the last 
few years (I]] . However physical implementations of large 
scale quantum networks are severely limited by decoher- 
ence g, the loss of coherence in the superposition state 
of the system. This is especially true of extended sys- 
tems that involve entanglement of macroscopic degrees 
of freedom, such as the center-of-mass motion of atoms 
in a trap || . 

We consider multiple computational levels in each 
atom to store information more densely. The radiative 
lifetime of Rydberg levels scales as n 3 , where n is the 
principal quantum number, and approaches a millisec- 
ond for n > 100. In most atomic quantum logic schemes, 
the bottleneck for the coherence time is the macroscopic 
entanglement of two or more atoms. The multilevel ap- 
proach minimizes this decoherence by reducing the num- 
ber of atoms in the quantum network, with each d-level 
atom representing log 2 d qubits of information [Q. In a 
d = 8 level scheme, for example, the number of atoms 
that need to be entangled is less by a factor of three. 

A great deal of research has been carried out in the past 
few years developing techniques for controlling multilevel 
Rydberg wave packets using shaped short pulses. It is at- 
tractive to exploit this work by developing an alternative 
wave packet basis for the quantum information stored in 
the atomic energy levels. Rydberg wave packets are co- 
herent superpositions of atomic energy eigenstates with 
large principal quantum numbers. We can use the rela- 
tive phases of the different eigenstates in the superposi- 
tion to store information. Using the energy-time Fourier 
kernel in unitary time evolution, we define an orthogonal, 
discrete wave packet basis that is related to the energy 
levels by a quantum Fourier transform. This transform 
plays a key role in quantum computing applications j^], 
and we show how it naturally permits complementary 
bases for representing the information stored in a multi- 
level quantum system. 

In the optical domain, a discrete Fourier transform 
connects the spectral and temporal modes of the ra- 
diation field in a cavity. Iaconis et al. M use this 



approach to characterize the Husimi distribution of a 
multimode quantum field in heterodyne detection. Tit- 
tel et al. Q| exploit the complementarity of energy and 
time in their scheme for quantum cryptography using 
entangled photons. In the atomic case, Muller and Noor- 
dam H introduce a Fourier relation between the energy 
quantum number and a continuous phase coordinate to 
study the photoexcitation of Rydberg states. We show 
how the quantum Fourier transform of the energy levels 
corresponds exactly to a discrete temporal basis of wave 
packets in cyclic motion about the nucleus. 

We consider radial wave packets, which are super- 
positions of Rydberg states with low angular momen- 
tum. These states are radially localized in probability 
and evolve periodically for short times, or in the classi- 
cal correspondence limit Q. The wave packet basis we 
introduce corresponds to discrete times in this classical 
evolution, when each wave packet state is centered at a 
different radius from the atomic nucleus. This basis is 
non-stationary in time, with wave function amplitudes 
evolving cyclically during a classical period. Due to the 
unequal spacing of the energy levels, the wave function 
undergoes dispersion, which leads to revivals and super 
revivals in the wave packet amplitudes. The classical 
motion and quantum revival structure of this free time 
evolution allows for the possibility of conveniently con- 
trolling the wave function at suitably chosen intervals of 
time. 

We propose using short optical pulses to couple coher- 
ently the Rydberg manifold to a low-lying energy level 
in the atom. A transform-limited pulse that is short 
compared to the classical Kepler period of the system 
interacts with all the energy levels but couples only to 
the wave packet state that is nearest the atomic core, a 
phenomenon familiar in the excitation and photoioniza- 
tion of Rydberg states (h]] . The classical electron has the 
largest momentum, and hence absorbs energy most effi- 
ciently, when near the core. During the pulse, the wave 
packet state nearest the nucleus and the ground state 
form an effectively two-level system, and undergo Rabi 
oscillations. 
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Coherent control of wave packets based on spatial lo- 
calization has been studied by Ahn et al. pj| , who show 
that half-cycle pulses can be used to extract binary phase 
information stored in a radial wave packet. The use of 
optical pulses shorter than a Kepler period for quantum 
control of atomic wave packets in the strong-field regime 
has been studied analytically by Araujo et al. jj^j . A 
weak field analysis by Noel and Stroud Jl3| showed that 
a uniform train of d pulses within one Kepler period can 
be used to create an arbitrary Rydberg state of d levels 
whose amplitudes are related by a discrete Fourier trans- 
form to the pulse amplitudes. We propose to use such 
pulses for selective processing of quantum information in 
the wave packet basis. 

Consider d Rydberg energy levels with angular mo- 
mentum I = 1 in a hydrogen or alkali atom as a quantum 
information basis, 



1^ = 171 + ^1,0), 

j = -d/2 + l,-d/2 + 2,. 



(1) 



,d/2, 



where n is the mean value for the principal quantum num- 
ber in the wave packet superposition, and the subscript 
v is used to denote a state in the energy basis. We have 
assumed above that the number of levels d is an even 
number, but the arguments below are easily extended to 
odd d. These levels can be excited simultaneously from 
the ground state by a strong linearly polarized laser pulse 
with sufficient bandwidth to overlap all the levels. The 
pulsed excitation creates a wave packet which evolves in 
time as 



(2) 



where tkJj is the energy of the jth Rydberg level and 
the amplitudes bj are determined by the spectrum of the 
exciting pulse. If we assume that the pulse spectrum, 
centered at ujq, is uniform over the d computational lev- 
els and falls off rapidly outside this manifold, we have 
a uniform amplitude distribution, bj = l/Vd, which is 
time independent in the absence of external fields. 

The radial wave packet was first studied by Parker and 
Stroud ||, who showed that it corresponds to a radially 
localized shell of probability distribution moving periodi- 
cally between the classical turning points for short times. 
Eventually, the uneven spacing in the energies of the lev- 
els in Eq. ([j]) leads to dispersions and revivals of the wave 
packet, corresponding to the time-dependent phases in 
Eq. (||) interfering destructively or constructively. The 
rates at which the phases evolve with respect to each 
other depend on the frequencies u>j relative to the av- 
erage frequency ujq, which can be expanded in a Taylor 
series in j = n — n, 
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FIG. 1. Illustration of conjugate bases: Energy eigenstates 
\j) v related by a Fourier transform to wave packet states \k) 
evenly distributed in time around a classical orbit. A radial 
wave packet basis can be thought of as an ensemble of such 
orbits with different orientations for the ellipses. 



where in atomic units, Tk = 2irn 3 is the classical Kepler 
period for round trip motion, T rcv /2! = 27m 4 /3 is when 
the wave packet revives first, T sr /3! = 7rn 5 /6 is when 
the wave packet super-revives first, etc. At n = 180 for 
example, T# — 0-89 ns, T rev ~ 106 ns, and T sr ~ 14 /zs. 
For large n, the different time scales are well separated in 
magnitude, and for a small number of levels d, the wave 
packet is more spread out in space and disperses more 
slowly in time. For d 2 <C n, we see from Eq. (||) that the 
revival dynamics of the wave packet can be neglected to 
good approximation during the first few Kepler periods, 
and Eq. (0) becomes 



-^^exp(-i27rit/T K ) 



\^K{t)) 



(4) 



which evolves periodically in time with period Tr-. 

The Fourier kernel in Eq. (^J) involves j and t, and 
can be thought of as underlying an uncertainty relation 
between energy and time, with more levels leading to 
greater spatio-temporal localization of the wave packet. 
Although the Rydberg energies are discrete, the time t 
for the wave packet varies continuously around the or- 
bit. However, if we consider d discrete times during any 
Kepler period, 



t k - k(T K /d), 
k = -d/2 + l,-d/2 + 2 > . 



(5) 



we can construct an alternate discrete wave packet basis 
\k) T for the system of energy levels in Eq. (lj) based on 
the quantum Fourier transform (QFT), 
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\k) T = I^Mffc)) 

= ^Y^e W (-i2njk/d)\j) u . (6) 

This is an orthogonal basis in time (denoted by the sub- 
script t), and the basis amplitudes evolve periodically 
with the classical period Tk, as illustrated in Fig. [|. 
Thus fc — and fc — d/2 correspond to wave packet 
states centered at the inner and outer turning points of 
the classical orbit respectively. The sign of k indicates 
the direction of propagation of the wave packet. We can 
also define the basis for odd d, except that in this case 
k = ±(d— 1)/2 correspond to wave packet states centered 
near (but not exactly at) the outer turning point, having 
nonzero and opposite momenta. The energy eigenstates 
can be recovered from \k) T by an inverse QFT, 

V d k 

In general, any excited state |^(f)) of the d Rydberg 
levels can be expanded in either basis, 

!*(*)> =!>(*) expH^tJIj),, (8) 
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= exp(~iuj Q t)Y,h(t)\k) T , (9) 
fc 

where we have removed the center frequency u>o from the 
wave packet amplitudes bk- A classical discrete Fourier 
transform (DFT) relates the two sets of amplitudes, 

bk(t) — —j= ^ bj(t) exp(-iojjot) exp (i2njk/d) . (10) 
V d 

j 

For free time evolution, bj(t) = bj(0), and Eq. ( |l0| ) can 
be written as 

bk(t)=J2~bk'(0)uk'-k(t), (11) 

fc' 

where Uk>- k (t) is the DFT of the time evolution kernel, 
u k '-k(t) = - ^exp [-iujjot - i2nj(k' - k)/d] . (12) 

3 

In the Kepler regime of Eq. (|J), the levels are equally 
spaced, Wjo — ^j/Tk, and Uk'-k becomes a delta func- 
tion in fc', resulting in a cyclic permutation of the wave 
packet amplitudes in Eq. (pi]), 

b k {nT K /d) = bk-niO), (13) 

where n is an integer and fc — n is taken modulo d. Equa- 
tion ( |l3| ) amounts to a SHIFT gate in the wave packet 
basis, the multivalued equivalent of the NOT gate for bi- 
nary logic. More complex unitary transforms occur for 
longer times, when the revival and superrevival dynamics 



become relevant. Thus free time evolution alone becomes 
useful for quantum information processing in the wave 
packet basis. This is a consequence of the nonstationary 
character of the basis and the unequally spaced energy 
levels. 

We now turn to the control of the multilevel atomic 
basis using laser fields. Universal multilevel transforms 
of the amplitudes bj can in principle be achieved us- 
ing narrow-band lasers to couple neighboring transitions 
among the energy levels Q. However this requires mul- 
tiple laser frequencies tuned to the d — 1 transitions. A 
simpler control scheme is possible for the wave packet am- 
plitudes as shown below. This involves short broad- 
band pulses and takes advantage of the cyclic evolution 
of the wave packet states described above. 

A laser pulse with a spectrum overlapping the d Ryd- 
berg levels couples most strongly to only that combina- 
tion of the levels with a phase relation between them that 
localizes the electron near the core. This corresponds to 
the wave packet state |0) r , for which all the phases are 
equal since fc = 0. Consider a pulsed field with a carrier 
frequency u> u> and pulse profile /(f), 

E(f) = eEof(t) exp(-iwf) + c.c, (14) 

where Eq and e are the amplitude and polarization of 
the field. Assume that the pulse couples the d Rydberg 
levels to the ground state \g) in the atom. In the energy 
level basis for the Rydberg states, the dipole interaction 
with the field yields the following coupled equations for 
time evolution in the rotating wave approximation, 

^^/(^^cxpHA.t)^, (15) 

3 

b j = lf(t)n j exp(iA j t)b g , (16) 

where Qj and Aj are the Rabi frequency and field detun- 
ing defined as 

Qj = 2{g\d\j) v ■ eE /h, (17) 
Aj = Uj - uj, (18) 

and we have taken fij to be real in Eqs. ( |l5| ) and Jl6|). 
Using the Fourier inverse of Eq. (0), we rewrite the equa- 
tion of motion for b g in terms of the wave packet ampli- 
tudes, 

b g = ^f(t)e X p(-iA t)J2&kbk, (19) 

fc 

where Ao is the mean detuning of the field, and is the 
DFT of the Rabi frequency fij, 

exp (— i2irjk/d) . (20) 

Vd - 

To find the equation of motion for the wave packet am- 
plitude bk, we take the time derivative of Eq. ( |l0| ) and 
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FIG. 2. Creation of a dark wave packet at n = 180 and 
d = 8. The spectral width (FWHM) of the n pulse is approxi- 
mately 1.3d/Tftr. The wave packet populations \bk\ 2 are shown 
before (t = -T K /d), during (t = 0) and after (t = T K /d) the 
pulse. 




FIG. 3. Wave packet control for quantum information pro- 
cessing. The energy level basis is \j) = \fi + j, 1,0). Broad- 
band 7r pulses de-excite two wave packet amplitudes in the 
conjugate wave packet basis. The ground levels g and e are 
coupled by a laser with frequency u ~ uo e . 



use Eq. ( |l6| ) for bj. Assuming that the pulse width is 
short compared to a Kepler period, we can ignore the free 
atomic evolution during the pulse that led to Eq. (JTT 
In the short pulse regime, we are left with 



b k = ^/(*) fi fe ex P()'Aoi)V 



(21) 



Equations ( |l9| ) and (^l|) describe Rabi oscillations be- 
tween the ground state and the wave packet states, where 
the Rabi frequencies f2fc are Fourier transforms of their 
energy level counterparts Qj. For Rydberg levels, flj 
varies slowly with j, 



tij ~ (n + j) 



-3/2 



(22) 



making only the dc component £!o significant for large n. 
For n = 180 and d = 8, we find 0±i ~ 0.01 Q . Thus to 
good approximation, we can ignore the Rabi frequencies 
f2fe/o in Eqs. ( |l9|) and (|2l]). These equations reduce to 



■f(t)& exp(-iA t)b , 



bo = ^ /(i)^oexp(zA i)6 g . 



(23) 
(24) 



These describe two-level Rabi oscillations between the 
ground state and the innermost wave packet state dur- 
ing the pulse, with a detuning Ao = loq — u> between 
the peak of the pulse spectrum and the average Rydberg 
frequency. The two-level approximation is valid only for 
short pulses. To affect only the k = wave packet state, 
we need the pulse width r p to be less than Tx/d, the time 
taken for this wave packet element to leave the atomic 
core. However, the bandwidth of the pulse cannot be 
much larger than the atomic frequencies spanning the 
d-level basis, to avoid coupling into levels outside this 
basis. Since the Rydberg frequency spacing is approxi- 
mately 1/Tk, this implies a bandwidth Auj/2n ~ d/Tx- 
To summarize, we require that 



Tk 
d 



2tt 
Aw' 



(25) 



or that the pulse be transform- limited. For n = 180 
and d = 8 levels, we use a Gaussian pulse of width 
t p = 0.5 ln2 T K /d ~ 110 ps (FWHM amplitude), cor- 
responding to a spectral width Aoj/2it — (8/2ir) d/TK 
for the electric field. Figure shows an exact numeri- 
cal solution for the wave packet amplitudes before, dur- 
ing and after the pulse. Assuming that the population 
is uniformly distributed initially, we find that a tt pulse 
transfers nearly all the population in the k = state to 
the ground state, creating a "dark" wave packet in the 
Rydberg basis that evolves cyclically in the Kepler evo- 
lution. About 5% of the population is lost from neigh- 
boring states k — ±1 during the pulse. Thus a coherent 
coupling of individual wave packet states with the ground 
state can be modeled as a two-level system to good ap- 
proximation. 

Since the wave packet basis is nonstationary in time, 
the Rydberg wave function undergoes dispersion and re- 
vivals as the energy level phases evolve in the Schrodingcr 
picture. This leads to aperiodicity in the free time evo- 
lution of the wave packet amplitudes bk, going beyond 
the cyclic regime of Eq. (p"3|). At n = 180 and d = 8 for 
example, we find that the population in a wave packet 
state decays by about 5% after a Kepler period, but re- 
covers to within 2-4% at a revival or super-revival. This 
scales as n~ 2 , making dispersion losses less significant 
for larger n. By appropriate timing of successive laser 
pulses, we can thus take advantage of the revival struc- 
ture of the atomic energy spectrum when addressing the 
wave packet amplitudes in the basis. 

An arbitrary unitary transform of the d wave packet 
states can be decomposed into a product of 0[d 2 ] two- 
dimensional transforms |l4"| , each of which can be imple- 
mented sequentially on two de-excited wave packet am- 
plitudes at a time, as shown in Fig. ||. Two low lying 
energy levels g and e temporarily store the wave packet 
amplitudes in the form b^l g) + 6/c'|e). Phase shifts or ro- 
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tations of this state are accomplished by coupling the two 
levels with a pulse of appropriate length, Rabi frequency 
and detuning. The transformed state is then restored 
to the Rydberg basis. Thus multilevel transforms in the 
atom are possible with a timed sequence of laser pulses, 
without the need for multiple lasers to address the differ- 
ent Rydberg transitions. Moreover the fidelity of the gate 
operations compares favorably with that in an equivalent 
system of entangled two- level atoms || . 

For efficiently scalable quantum computing, we re- 
quire conditional transforms on two entangled multilevel 
atoms. The method used in Ref. || is based on the lin- 
ear ion trap scheme for binary quantum logic pq |. This 
scheme provides a means to evolve one two-level atom in 
the trap array conditional on the state of another two- 
level atom. Since the wave packet approach allows each 
basis amplitude in the atom to be addressed sequentially 
through low-lying energy levels, conditional transforms 
on other atoms can be carried out using the same meth- 
ods as in the binary scheme. That is, conditional on a 
given wave packet state (vis a vis energy level) in the con- 
trol atom, a d-dimensional transform is performed on the 
wave packet basis of the target atom using the methods 
outlined above. 

The time taken to couple two atoms in this approach 
can be on the order of the atomic Kepler period, which 
permits a fast non-adiabatic coupling scheme. Also, the 
mean energy and number of Rydberg levels used depend 
ultimately on the ionization threshold imposed by exter- 
nal fields in the environment. We anticipate that neutral 
atom schemes that allow fast coupling times, as have been 
suggested recently |16 , will provide a useful setting for 
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